Chemical Equilibrium in Collisions of Small Systems 
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The system-size dependence of particle production in heavy-ion collisions at the top SPS energy 
is analyzed in terms of the statistical model. A systematic comparison is made of two suppres- 
sion mechanisms that quantify strange particle yields in ultra-relativistic heavy-ion collisions: the 
canonical model with strangeness correlation radius determined from the data and the model formu- 
lated in the canonical ensemble using chemical off-equilibrium strangeness suppression factor. The 
system-size dependence of the correlation radius and the thermal parameters are obtained for p-p, 
C-C, Si-Si and Pb-Pb collisions at ^/sNN = 17.3 ^GeV. It is shown that on the basis of a consistent 
set of data there is no clear difference between the two suppression patterns. In the present study 
the strangeness correlation radius was found to exhibit a rather weak dependence on the system 
size. 
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I. INTRODUCTION 

Experiments with ultra-relativistic nucleus- nucleus 
collisions provide a unique opportunity to study the prop- 
erties of strongly interacting matter under the extreme 
condition of high energy densities. Hadronic multiplic- 
ities, and their spectra in particular, carry information 
about the nature of the medium from which they origi- 
nated. The statistical model (SM) has been recognized 
as a powerful approach to describe particle production 
yields in heavy-ion collisions [l], These models as- 
sume that particle creation occurs at chemical freeze-out 
with the collision fireball being in thermal and chemical 
equilibrium. 

In the limit of high temperature and/or large system 
size, the grand-canonical (GC) treatment of strangeness 
conservation is adequate in the SM. There, strangeness 
conservation is implemented on the average and is con- 
trolled by the chemical potential [l], However, if 
the number of strange particles in a collision fireball is 
small, either due to low temperature (e.g. in A-A col- 
lisions at incident energies of a few ^GeV 0) or due 
to small system size (e.g. in a fireball created in p- 
p or light-ion collisions), then strangeness conservation 
has to be implemented exactly in the canonical (C) en- 
semble [H, S H, H S H, Hi- Exact strangeness conser- 
vation leads to the suppression of the strange particle's 
phase-space and is usually referred to as canonical sup- 
pression [1, H, However, canonical suppression with 
the assumption of strangeness chemical equilibrium in 
the whole fireball volume was found to be insufficient 
to reproduce the observed yields in nucleus- nucleus col- 
lisions at the SPS |1, i, [lo|, [H [l2|. In this paper we 
consider a method to account for the suppression beyond 



the one expected in the canonical ensemble by assuming 
that exact strangeness conservation holds only in a small 
sub- volume Vc of the system [H, 0: 0| ■ The concept of 
such a sub-volume, or a strangeness correlation volume, 
has been used in earher studies 0, [H, 0, [11] ■ Here, we 
present a systematic analysis of p-p and central C-C, Si- 
Si and Pb-Pb collisions at the top SPS energy from the 
NA49 collaboration [H, [TJ [H, [ll, [H M, 
We selected a consistent set of data (discussed below) 
for different collision systems and focus only on central 
nucleus-nucleus collisions. We test the validity of the 
above approach and study the relation between Vc and 
the system size, and also the chemical freeze-out volume. 
The results of the canonical model with cluster forma- 
tion will be compared to those obtained in the canonical 
model with strangeness suppression through chemical off- 
equilibrium factor 7s. 

The paper is organized as follows: A short outline of 
the statistical model and its canonical formulation with 
cluster formation is given in Section [Til In Section [ml 
the data analysis is presented followed by a comparison 
of the models with the data in Section IIVI The paper 
concludes in Section [Vl with a discussion. 



II. THE MODEL DESCRIPTION 

A. Grand-canonical ensemble 

In the statistical model the thermal fireball is char- 
acterized by its volume V, temperature T and charge 
chemical potentials fl which are assumed to be uniform 
over the whole volume. In the GC ensemble, these pa- 
rameters determine the partition function Z{T,V, fl). In 
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FIG. 1: Different particle ratios as a function of the radius R 
of a spherical volume. The temperature T = 170 MeV and 
the baryon chemical potential /is = 255 MeV were chosen 
according to the thermal conditions at top SPS energy. All 
ratios are normalized to their grand-canonical values. 



the hadronic fireball of non-interacting particles and res- 
onances, h\Z is the sum of the contributions of all i- 
particle species of energy Ei and spin-isospin degeneracy 
gi by, 



\nZ(T,V,ll)=Y,Zl{T,il), 



(1) 



to hadron i with the branching fraction Tj^i have to be 
calculated as. 



n 



decay 



(4) 



Consequently, the final yield Ni of particle species i is the 
sum of the thermally produced particles and the decay 
products of resonances. 



N,, 



{n.^+n^ ')V. 



(5) 



From Eqs. ([3][5]) it is clear that in the GC ensemble the 
particle yields are determined by the volume of the fire- 
ball, its temperature and the baryon chemical potential. 



B. Canonical ensemble 

In general, if the number of particles carrying quantum 
numbers related to a conservation law is small, then the 
grand-canonical description no longer holds. In such a 
case conservation of charges has to be implemented ex- 
actly in the canonical ensemble. Here, we refer only to 
strangeness conservation and consider charge and baryon 
number conservation to be fulfilled on the average in the 
GC ensemble because the number of charged particles 
and baryons is much larger than that of strange particles. 
The density of strange particle i carrying strangeness s 
can be obtained in the canonical ensemble from, 



with. 



^'=^^Jo ^^dWn[l±exp(^^ 



(2) 



where ql — {Bi,Si,Qi) are the quantum numbers of 
particle i and fl — [fis, (J-St I^q) a-re the chemical po- 
tentials related to the conservation of baryon number, 
strangeness and electric charge, respectively. The upper 
(lower) signs refer to fermions (bosons). 

In a fireball created in heavy-ion collisions only the 
baryon chemical potential fiB is an independent pa- 
rameter while the other two, the charge hq and the 
strangeness /is chemical potentials, are constrained by 
the initial conditions from the electric charge of the in- 
coming nuclei and the condition of strangeness neutrality. 

The partition function ([1]) contains all information to 
obtain the number density rii of thermal particle species 
i. Introducing the particle's specific chemical potential 
/ij, one gets. 



n,(T, fl) 



1 djTlnZ) 



(3) 



The thermally produced resonances that decay into 
species i contribute to the measured yield. Therefore, 
the contributions from all heavier particles j that decay 
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where Z^'^q is the canonical partition function 



(6) 
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and Z} is the one-particle partition function ([2]) calcu- 
lated for — Q in the Boltzmann approximation. The 
arguments of the Bessel functions Ig {x) and the parame- 
ters ai are introduced as. 



2V./SIS- 



(8) 



where Ss is the sum of all [fis = 0) from ^ for particle 
species k carrying strangeness s. 

In the limit where x„ < 1 (for n = 1, 2 and 3) the 
density of strange particles ^ carrying strangeness s is 
well approximated by 



' Io{xi) ' 



(9) 
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From this equation it is clear that in the C ensemble the 
strange particle density depends explicitly on the vol- 
ume V through the arguments of the Bessel functions. 
In addition, as seen in Eq. the strange particle den- 
sities are suppressed in the C ensemble due to exact 
strangeness conservation constraints. 

The canonical suppression relative to the GC results 
for different particle ratios is illustrated in Fig.[l]for fixed 
T and /ig and for different system sizes parameterized 
by the radius R, assuming spherical geometry of the vol- 
ume V. The left-hand panel illustrates that the canonical 
suppression increases with the strangeness content of the 
particle For i? > 5 fm the canonical suppression is 
already so small that the C and GC descriptions of par- 
ticle ratios are equivalent. The right-hand panel shows 
the particle ratios with a difference in their strangeness 
quantum numbers, AS = 1. Interestingly, the canonical 
suppression of these ratios is similar but not identical. 
Thus, the C-suppression in these ratios depends on the 
strangeness content of both hadrons and not only on the 
difference. 



C. Strangeness suppression mechanisms 

In the application of the statistical model to particle 
production in heavy-ion and particularly in elementary 
particle collisions it was found that canonical suppression 
alone is not sufficient to quantify the observed strange 
particle yields |6|. Consequently, additional suppression 
mechanisms were proposed to account for deviations from 
experimental data. Here we present two methods that 
lead to additional suppression of strange particlcphase- 
space going beyond the normal canonical effect [J 0, 0] • 

The suppression of strangeness has been parameterized 
by a factor, 7s, that is introduced to suppress hadrons 
composed of strange and/or anti-strange quarks [lQ|. In 
this description the strange particle density ([3]) composed 
of s strange quarks/antiquarks is modified in the GC en- 
semble by, 



7s "-i, 



(10) 



where 75 is an additional parameter of the model. This 
chemical off-equilibrium factor also modifies the calcula- 
tion of the canonical suppression outlined in Eqs. ([Ml]), 
the single-particle partition function is replaced by 

In order to account for an additional strangeness sup- 
pression in terms of the canonical model we use the con- 
cept of strangeness correlation in clusters of sub- volume 
Vc < 1^ [H . Consequently, there are two volume param- 
eters in the model; the overall volume of the system V, 
which determines the particle yields at fixed density and 
the strangeness correlation (cluster) volume Vc, which 
enters through the canonical suppression factor and re- 
duces the densities of strange particles. Hence, this re- 
sults in replacing the volume V by Vc in Eq. As- 
suming spherical geometry, the volume Vc is parameter- 
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FIG. 2: Midrapidity (circles) and Att (squares) data on the 
A/A and K~/7r~ ratios in p-p and central C-C, Si-Si and Pb- 
Pb collisions normalized to the Pb-Pb measurement [l6l . [itI . 



ized by the radius Rc which serves as a free parameter 
and defines the range of local strangeness equilibrium. A 
particle with strangeness quantum number s can appear 
anywhere in the volume V. However, it has to be ac- 
companied within the sub-volume Vc by other particles 
carrying strangeness — s to conserve strangeness exactly. 

In the application of the SM to particle production and 
their system-size dependence we consider three different 
formulations of the model: 

(a) equilibrium model in the canonical ensemble with 

strangeness correlation volume Vc — V, 

(b) non-equilibrium canonical model with strangeness 

undersaturation parameterized by the factor 75, 

(c) canonical model that accounts for strangeness corre- 

lation in the sub- volume Vc < V that is quantified 

by the cluster radius Re- 
in the following we compare the above models with ex- 
perimental data at the top SPS energy for different col- 
liding systems. We focus on the system-size dependence 
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TABLE I; Particle yields (47r integrated) in minimum bias p- 
p collisions at = 17.3 GeV from Ref. [13, \^ and fit 

results from model (c). 



Particle 


Ti-rt 
P F 

Data 


Fit 




1.23±0.26 


1.23 


K+/K- 


1.61±0.61 


1.65 


A/A 


0.118±0.013 


0.118 


K-/-JT- 


0.062±0.021 


0.062 


A/K+ 


0.450±0.138 


0.441 


A/K- 


0.085±0.027 


0.086 




2.57±0.36 


2.57 



of the thermal parameters with particular emphasis on 
the change in the strangeness correlation radius Re- 



III. DATA SETS 

To verify the statistical nature of particle production 
in heavy-ion collisions for different colliding systems we 
consider the experimental data at the top SPS energy 
with y/sNN = 17.3 AGeV, froin p-p and central C-C, Si- 
Si and Pb-Pb collisions [11, M, M, M,M,M,M,m- 
In view of the fact that only a limited number of 
hadrons were experimentally analyzed in the smaller 
systems, we restrict our study to a consistent set of 
data [11, [13, HiH, Hi. Therefore, e.g. the multi-strange 
particles are not included as they are not available for 
all data sets. Both midrapidity densities and integrated 
yields are studied with the exception of the p-p collisions 
where the midrapidity data are not available. A compi- 
lation of p-p and heavy-ion data used in our analysis is 
summarized in Tables HI and [Til 

The system-size dependence of ratios of particle yields 
is illustrated in Fig. [2] for the K^/tt^ and A/A. The 
K^/tt^ ratio shows rather moderate dependence on the 
system size. The midrapidity data on the K~/7r~ ra- 
tio are almost saturated at the value expected for cen- 
tral Pb-Pb collisions already for A^part > 20, whereas the 
rapidity-integrated ratio shows such a property only for 



A'part > 30. The variation of the A/ A ratio with the 
system size is stronger than that observed in the K~ /tt^ 
ratio. At mid-rapidity the A/A ratio exceeds its value in 
the full phase-space by almost a factor of two. 

The yield of (/)-mesons per pion in different colliding 
systems, normalized to its value in the most central Pb- 
Pb collisions, is shown in Fig. |3l The yield of (/)-mesons 
is seen to decrease stronger with decreasing system size 
than the pion yields. This indicates that the 4> with 
hidden strangeness does not behave like a strangeness- 
neutral object. For comparison, also shown in Fig. |3l are 
the results for the system-size dependence of the normal- 
ized K/tt and (K/tt)^ ratios. From the comparison with 
kaons its is not clear whether the A/tt ratio favors the 
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FIG. 3: The K/tt, {K/nf and (jif-jT ratios of Att integrated 
yields normalized to their values in the central Pb-Pb colli- 
sions at the SPS. The p-p and central C-C, Si-Si and Pb-Pb 
data are from NA49 Collaboration at ^/snn ~ 17.3 AGeV 
[H, [13, [H, miH, m, m. The open symbols refer to non- 
central Pb-Pb collisions. 



single- or double-strange particle properties. Instead, it 
seems to behave as a strange particle with an effective 
quantum number between one and two. 



IV. MODEL COMPARISON WITH DATA 

In the SM the particle production yields are fully deter- 
mined by the thermal conditions in the collision fireball. 
In all models (a), (b) and (c) there is a common set of 
parameters {T,iib,V) that quantifies the thermal parti- 
cle yields. In the SM models (b) and (c), that account 
for an additional strangeness suppression, there are ad- 
ditional free parameters: 75 and Rc, respectively. For 
our study, we used the code THERMUS [25 1. 

Before considering a detailed quantitative analysis of 
particle production for different systems within the SM 
we first discuss general trends expected within the models 
based on the data shown in Figs. [2] and [3l 



A. General trends 

In terms of the models (a), (b) and (c) the anti- 
strange to strange particle ratios are independent of the 
fireball volume. In the C ensemble an explicit depen- 
dence of Us on the volume appears in the same way for 
particles and anti-particles, thus it is cancelled in their 
ratios. Consequently, the strong variation of the A/A ra- 
tio with the system size seen in Fig. ^ could be related 
to the variation in T and/or fj,B with the system size. It 
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TABLE II: 47r-integrated particle yields (upper part) and midrapidity densities (lower part) in central C-C, Si-Si and Pb-Pb 
collisions at ^/snn ~ 17.3 GeV and fit results from model (c). The experimental results were taken from Ref. [IB, El for C-C 
and Si-Si and from Ref. for Pb-Pb. 





C-C 




Si-Si 




Pb-Pb 




J. CLl. LlVjlt. 


Data 


Fit 


Data 


Fit 


Data 


Fit 






1.01 


n Qs + 1 1 

\J.iJKJ I \J . A. A. 


0.98 


97 + 10 


0.99 


/K~ 


1 70 ± 30 


1.71 


1 68 ± 26 


1.70 


1 98 ± 27 

A- m 1 V7 « ^ 1 


1.77 


A/ A 


\J . _1_ \J.\JOiJ 


0.134 


1 97 + n 033 


0.127 


n 083 + 093 


0.080 


K~ /tt~ 


067 ± Oil 


0.067 


077 ± 010 

\J m\J t 1 — 1 — \J m\J A\J 


0.077 


081 ± 010 

KJ A. — 1 — \J t\J A-yj 


0.082 


\/K^ 




0.52 


52 + 1 3 


0.51 


n 44 + 1 

yjw^^ 1 \J , AVJ 


0.51 


A/K~ 


119 ± 031 

\J M A A tj — 1 — \J M\Jt^ A 


0.119 


111 ± 036 

\j , A A A — 1 — \J Myj^jyj 


0.111 


072 ± 015 


0.072 




22.2 ± 1.9 


22.2 


57.6 ± 4.6 


57.6 


639 ± 48 


639 




0.98 ± 0.21 


1.07 


0.99 ± 0.19 


1.03 


0.97 ± 0.08 


0.99 


K+/K- 


1.60 ± 0.28 


1.30 


1.60 ± 0.26 


1.46 


1.76 ± 0.15 


1.62 


A/A 


0.269 ± 0.085 


0.276 


0.182 ± 0.070 


0.179 


0.140 ± 0.044 


0.122 


K~/tt~ 


0.093 ± 0.018 


0.094 


0.100 ± 0.018 


0.101 


0.096 ± 0.008 


0.097 


A/K+ 


0.280 ± 0.063 


0.312 


0.366 ± 0.081 


0.396 


0.361 ± 0.082 


0.459 


A/R- 


0.121 ± 0.034 


0.112 


0.107 ± 0.037 


0.104 


0.089 ± 0.022 


0.091 




5.70 ± 0.85 


5.70 


15.0 ± 2.0 


15.0 


175.4 ± 9.7 


175.4 



is quite natural to assume that the temperature reached 
in the coUision scales approximately with the collision 
energy. Therefore, at fixed ^/snn the temperature is ex- 
pected to show a weak variation with the system size. On 
the other hand, the chemical potential can be viewed as 
the measure of baryonic stopping in the collision. Thus 
it should decrease with ^/snn and increase with system 
size. With the above generic properties of thermal pa- 
rameters, the strong variation of A/ A with the system 
size seen in Fig. [2] could be attributed to a change in /i^ 
at T ~ const., since in the SM A/A ~ exp {-2fiB/T). 

In the SM the K~ yield is not explicitly dependent 
on fiB- However, due to strangeness neutrality, /is = 
/^s(/^s), implying a weak influence of fiB on the kaon 
yield. Consequently, the change in the K/tt ratio with 
the system size is mainly to be expected as a consequence 
of canonical suppression and/or chemical off-equilibrium 
effects. 

In the equilibrium statistical models (a) and (c) the (/)- 
meson, being a strangeness neutral particle, should show 
similar system-size dependence as any other non-strange 
meson. Thus, the strong variation of the (jy/n ratio with 
centrality seen in Fig. [3] is inconsistent with the SM (a) 
and (c). The system-size dependence of </()- mesons is 
only expected to be qualitatively different from the other 
mesons in the model (b) that accounts for strangeness 
suppression by including chemical off-equilibrium param- 
eter 75. In this model the 0-meson, being made of a 
strange/antistrange quark pair, is suppressed by the fac- 
tor 7I, thus behaving as a strangeness- two particle. How- 
ever, from the discussion in the last section it is clear that, 
on the basis of data, it is not possible to fix the effective 
strangeness content of the i^-meson. Consequently, to 
avoid ambiguities, we have excluded the hidden-strange 
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FIG. 4: Midrapidity 7r~-density and midrapidity particle ra- 
tios from C-C collisions at y^ijvjv = 17.3 AGeV from NA49 
0, together with different model predictions: (a), (b) 
and (c) introduced in the text. The lower panel shows the 
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particles from the SM analysis. 
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FIG. 5: Integrated tt yield and integrated particle ratios 
from Pb-Pb collisions at ^/i]^ = 17.3 AGeY together 
with the model fits as in Fig. |4l The lower panel shows the 
deviation of the model fits to data. 



B. System-size dependence of thermal parameters 

The bulk properties of the system-size dependence of 
different particle yields in A-A collisions at the top SPS 
energy show qualitative agreement, except for the 0- 
meson yields, with the expectations of the statistical 
models. To verify the validity of the models requires 
detailed quantitative studies. In the following, we ana- 
lyze the SPS data from p-p, C-C, Si-Si and Pb-Pb col- 
lisions in terms of three different concepts that were in- 
troduced as: (a), (b) and (c) in the previous section. 
The investigated yields and ratios (listed in Tables [J and 
HI)) together with model predictions are shown in Figs. 2] 
and [5] for C-C and Pb-Pb collisions respectively. Par- 
ticle ratios from the fits with the best ^re displayed 
for all of the models. The model (a) with a completely 
equilibrated strangeness abundance exhibits large devia- 
tions from data. The normalized differences between the 
model results and the data are shown in the lower panels 
of these figures. For the small systems, model (a) fea- 
tures very large values of ~ 2 per degree of freedom. 
Thus, the equilibrium SM with strangeness conservation 
in the C ensemble is not favored by the data. In contrast, 
both models that allow for an extra suppression of the 
strange-particle phase-space with 7s or with agree 
quite well with the measurements and yield comparably 
good descriptions of data. The latter can be correlated 
to properties of the hot and dense medium and system- 
atic studies might provide access to the production and 
hadronization mechanism of strange particles. 

The system-size dependence of thermal parameters ob- 
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FIG. 6: Chemical freeze-out temperature T and baryon chem- 
ical potential fiB from fits to midrapidity densities (circles) 
and integrated yields (squares) from p-p and central C-C, Si- 
Si and Pb-Pb collisions obtained in the canonical model (c) 
introduced in the text. The open circles refer to the baryon 
chemical potential at midrapidity extracted from the A/ A ra- 
tio keeping T as for C-C collisions. 



tained in these studies within the models (a) and (b) 
agree with previously published results and thus will not 
be discussed here. Instead, in Figs. [6l [7] and [8] we concen- 
trate on the system-size dependence of thermal parame- 
ters characterizing the canonical model (c) with cluster 
formation. The corresponding values are summarized in 

Table nn 

The freeze-out temperature T extracted from midra- 
pidity data is barely dependent on the system size, as 
shown in Fig. [5) This is in agreement with the expecta- 
tion that the chemical freeze-out temperature is predom- 
inantly established by the collision energy. Within sta- 
tistical errors the An data are also consistent with having 
the same chemical freeze-out temperature for all colliding 
systems from p-p up to Pb-Pb. 

No significant change in the baryon chemical potential 
is seen in Fig. [Hlfor An data. The variation of A/A seen 
in the Att data in Fig. [2] is a combined effect of minor 
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TABLE III: Statistical model parameters, extracted from the comparison of model (c) with experimental data from full phase- 
space (47r integrated) (upper part) and midrapidity densities (lower part) in minimum bias p-p, central C-C, Si-Si and Pb-Pb 
collisions at ^sjvjv ~ 17.3 GeV. 



Reaction 


r(MeV) 


Ats(MeV) 


_R(fm) 


i?c (fm) 


p-p 


174.7±21.2 


234±51 


1.31±0.39 


0.83±0.38 


C-C 


191.0±16.5 


275±53 


2.09±0.49 


0.74±0.20 


Si-Si 


181.9±14.6 


264±42 


3.19±0.63 


0.92±0.20 


Pb-Pb 


164.6± 7.5 


274±34 


8.9±1.1 


1.26±0.19 


C-C 


163.4±11.1 


128±29 


1.99±0.32 


1.22±0.30 


Si-Si 


166.3±11.8 


183±36 


2.58±0.46 


1.29±0.33 


Pb-Pb 


166.4± 7.8 


231±28 


5.71±0.67 


1.32±0.20 



changes in the freeze-out temperature and the chemical 
potential in different systems. At midrapidity, the chem- 
ical potential shows a decrease by almost 100 MeV from 
central Pb-Pb to C-C collisions at constant temperature. 
Consequently, there is a steep increase of the A/ A ratio, 
as seen in Fig. O 

From the above results it is clear that at y = the 
collision fireball created in A-A collisions appears at con- 
stant temperature but with varying /is that decreases 
with decreasing iVpart- For the rapidity integrated data 
such systematics cannot be concluded. 

In the canonical model (c) there are two volume scales 
that characterize the system: the fireball radius at freeze- 
out R and the strangeness correlation radius Rc- The 
fireball radius R is determined by the particle multiplic- 
ity e.g. the pion yield. Thus, a smaller pion density at 
midrapidity causes a smaller freeze-out radius than for 
the integrated yields. 

The strangeness correlation volume Rc is extracted 
from the system-size dependence of the strange to non- 
strange particle ratios. The resulting system-size depen- 
dence of Rc is shown in Fig. \7\ An increase of single- 
strange to non-strange particle ratios from p-p to Pb-Pb 
reactions is reflected in a variation of the strangeness cor- 
relation volume with the system size. The cluster radius 
varies between 0.7 fm and 1.3 fm. A larger radius in the 
fits of midrapidity data is correlated with the experimen- 
tal observation of larger /t:^ ratios at midrapidity as 
compared to Att yields. 



V. DISCUSSION 

We have already discussed in Fig. [1] that for 
strangeness-one particles the canonical suppression for 
R > (2 — 3) fm shows a very flat change with R and 
is hardly distinguishable from the asymptotic GC-value. 
Consequently, an exact determination of Rc based only 
on single-strange particles is very difficult and requires 
high precision data. For i? > 3 fm an increase of Rc by 
a large factor will not influence the particle yields sig- 
nificantly. Therefore, the results shown in Figs. [7] and [8] 
for TVpart > 30 or i? > 5 fm should be considered as a 
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FIG. 7: Cluster radius Rc as a function of the system size 
in A-A collisions from fits of the model (c) to midrapidity 
densities (circles) and integrated yields (squares) from p-p 
and central C-C, Si-Si and Pb-Pb collisions. 



lower limit. The actual change of Rc with system size 
and the chemical freeze-out radius could be stronger than 
that shown in Figs. [7|and [51 This happens, for exam- 
ple, at the SIS energy where the correlation volume Vc 
was shown to scale with the number of projectile partic- 
ipants 3 . A more precise determination of the relation 
between the volume at chemical freeze-out, system size 
and strangeness correlation at the SPS would require the 
analysis of the system size and/or centrality dependence 
of multi-strange particles in A-A collisions that show a 
stronger sensitivity to the strangeness correlation volume 

Relating the cluster to the fireball radii Rc — Rc{R), 
it becomes clear that Rc has a significantly weaker 
system-size dependence than R, as seen in Fig. [S] In 
p-p collisions the strangeness correlation length Rc is al- 
most as large as the fireball radius. In contrast, with 
increasing system size, the cluster volume can be even 
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FIG. 8: Cluster radius Rc as a function of the fireball radius 
R from fits of the model (c) to midrapidity densities (circles) 
and integrated yields (squares) from p-p and central C-C, 
Si-Si and Pb-Pb collisions. The broken lines show various 
choices of the Rc = cR°' relation. The full line indicates the 
proportionality Rc = R- 

a few times smaller than the entire fireball. The depen- 
dence of Rc on R shown in Fig. [5] can be parameterized 
as Rc = c ■ R" with a < 1 and c ~ 0.7 fm. From the 
consistent data set used in these studies the value of a ap- 
pears in the range 1/3 < a < 1/2. From this analysis we 
could conclude that the transition from suppressed (in p- 
p) to almost saturated (in Pb-Pb) single-strange particle 
production happens in a narrow range of the correlation 
length. This could indicate that Rc — (1-2) fm is the 



length scale of strangeness correlation. This suggestion 
gets implicit support from the percolation model calcu- 
lations [28l |. 

The appearance of a particle correlation volume in p- 
p collisions, introduced here for strangeness, comes also 
from other observables. The mass shift of the /5°-meson 
in the 7r"'"7r~ decay channel observed in p-p collisions (29l | 
could not be substantially accounted for by the phase- 
space distortion of a Breit-Wigner line shape. It was 
argued [soj that re-scattering and re-interaction of pions 
contribute an additional term to the p'^ spectral function. 
The best agreement with the measurement is achieved if 
a range of tt+tt" scattering of 0.73 fm is chosen. This 
value is in agreement with the strangeness correlation 
length Rc obtained in this study from the statistical 
model analysis of p-p data at the SPS. 

In summary, two modified versions of the canonical 
statistical model have been investigated. The experimen- 
tally observed strong suppression of the strange-particle 
phase-space is retraced to local strangeness conservation 
within small correlated clusters in the fireball. These 
sub- volumes, in consequence, cause strong canonical sup- 
pression and allow one to successfully reproduce the ex- 
perimental data. Furthermore, the cluster size is found 
to be only weakly dependent on system size. In all data 
under study, it is of the order of 1 to 2 fm. 
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